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An analysis of the structural properties of three-dimensional Coulomb clusters confined in a 
spherical parabolic trap is presented. Based on extensive high accuracy computer simulations the 
shell configurations and energies for particle numbers in the range 60 < N < 160 are reported. 
Further, the intrashell symmetry and the lowest metastable configurations are analyzed for small 
clusters and a novel type of excited state which does not involve a change of shell configuration is 
identified. 

PACS numbers: 52.27Gr,82.70.Dd 



Spatially confined mesoscopic charged particle systems have a number of unique properties not observed in con- 
ventional quasineutral macroscopic plasmas of electrons and ions in discharges or solids, electrons and holes in highly 
excited semiconductors and so on. With the help of confinement potentials it has now become routine to trap, for 
long periods of time, plasmas of a single charge (nonneutral plasmas), e.g. electrons and ions and even positrons in 
Paul and Penning traps 0,0) for an overview see Q, or colloidal (dusty) plasmas in discharge chambers, e.g. 0. By 
varying the confinement strength researchers have achieved liquid behavior and even Coulomb crystallization of ions 
|{| and dust particles [(J • These strong correlation phenomena are of exceptional current interest in a large variety 
of fields ranging from astrophysics (interior of Giant planets), high-power laser compressed laboratory plasmas, to 
condensed matter and quantum dots @ etc. Coulomb (Wigner) crystals are expected to exist in many White Dwarf 
stars. 

A particular property of trapped mesoscopic (N < 1, 000) clusters in spherical traps is the occurence of concentric 
shells with characteristic occupation numbers, shell closures and unusual stable "magic" configurations. Due to their 
close similarity to nuclei, metal clusters or atoms, these mesoscopic systems have been called "artificial atoms". A 
number of papers has been devoted to the exploration of the energetically lowest shell configuration (ground state) 
and metastable ("excited") states of two-dimensional (2D) artificial atoms, e.g. IToL fTT| and references therein. 

Recently, Arp et al. succeeded in the first experimental creation of spherically symmetric 3D clusters, so-called 
"dust Coulomb balls" ^2 an d presented an almost perfect crystal of N = 190 dust particles. This raises the question 
about theoretical configurations of mesoscopic 3D Coulomb balls which is the subject of this paper. It is natural 
to start with an analysis of the ground state and lowest metatstable states, referring finite temperature and melting 
properties, e.g. 01, to a subsequent study. 

The theoretical analysis of 3D artificial atoms is much more involved than in 2D and has so far mostly been 
restricted to small cluster sizes with often conflicting results, e.g. [l4l jlBl fl6| and references therein. Rafac et al. [l5| . 
correcting earlier results, identified the first shell closure at N — 12 (the 13th particle is the first to occupy a second 
shell) and presented detailed data, including ground state energies for N < 27, but they missed the onset of the third 
shell, as did Hasse et al. 0|. Tsuruta et al. extended the table to A" = 59 0|. The most extensive data, for up 
to a few thousand particles, has been presented by Hasse et al. and is a valuable reference for theoretical and 
experimental groups. However, their tables report excited states rather than the ground states for N — 28 — 31, 44, 54 
and practically for all N > 63 (except for N = 66). The reason for the computational difficulties is the existence of 
a large number of excited (metastable) states which are energetically close to the ground state; with increasing N 
this number grows exponentially whereas the energy difference rapidly vanishes. This has to be accounted for by the 
computational strategy and choice of accuracy, see below. 

Model: We consider N classicle particles with equal charge q and mass m interacting via the Coulomb force and 
being confined in a 3D isotropic harmonic trap with frequency u> with the hamiltonian 



Despite its simplicity, model Q captures the basic properties of a multitude of classical systems, in particular of dust 
Coulomb balls and ions in spherical traps. Below we will use dimensionless lengths and energies by introducing the 
units r = (q 2 /2-nemuj 2 ) 1 ^ and E — (mw 2 q 4 /327r 2 e 2 ) 1 / 3 , respectively. 

To find the ground and metastable states, we used classical molecular dynamics (MD) together with a suitable 
"simulated annealing" method. Starting with a random initial configuration of N particles, the system is cooled 
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FIG. 1: Binding energy A2 (right axis) and mean Voronoi symmetry parameter (MVSP, left axis) for the two outermost cluster 
shells, a: N < 80, b: 80 < N < 160 (Color online). 

continuously until all momenta are zero and the particles settle in minima of the potential energy surface. Depending 
on the particle number, the cooling down process was repeated between a several hundred and a several thousand 
times until every of the computed low energy states was found more than a given number of times (typically 10 . . . 100) 
assuring a high probability (though no general proof) that also the ground state has been found. Crucial for a high 
search efficiency is the use of an optimized MD time step (it has to be chosen not too small to avoid trapping in local 
potential minima). The results are shown in tables HI and HT1 

Consider first the ground state shell configurations beyond the previously reported results see Tab. [H] 

Closure of the second shell is observed twice: for N — 57 16] and N — 60. Further, we have found the closure of 
the third shell to occur at N = 154, all larger clusters have at least four shells (in the ground state). The "nobel 
gas like" closed shell clusters are particularly stable, but a few others also have a comparatively high binding energy 
(addition energy change) A 2 {N) = E{N + 1) + E(N - 1) - 2E(N). Tsuruta et al. O found the stable clusters 
N = 4,6,10,12,19,32,38,56. For larger clusters the binding energy decreases, and the relative differences rapidly 
decrease. We found the next particularly stable ones to be N = 81, 94, 103, 116. The results are shown in Fig.^ The 
relative stability of these clusters is linked to a particularly symmetric particle arrangement within the shells which 
will be analyzed below. 

The existence of the shell structure is a marked difference from macroscopic Coulomb systems (N — > 00) and is, 
of course, caused by the spherical confinement potential. With increasing N the structure of a macroscopic system 
emerges gradually, see also Ref. 01 . This can be seen from the relative widths a m = a m /r m of the m— th shell 
(cr m denotes the variance of the shell radius r m ). For example, for N = 149, (starting from the outermost shell) 
o 1 = 0.0089, ct 2 = 0.035 and a 3 = 0.032, whereas for N = 160 we obtain o x = 0.0091, a 2 = 0.033 and a 3 = 0.0038. In 
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TABLE I: Shell configurations, energy per particle for the lowest lying states (for the excited states the energy difference with 
respect to the ground state is given) , mean radius of outer shell n , symmetry parameter Gm and number of Voronoi M-polygons 
TV(M) in brackets. For TV = 4, TV(3) = 4, and for TV = 5, N{3) = 2, TV (4) = 3. 



TV 


Config. 


E/N 


ri 


G 4 [TV(4)] 


G 5 [TV(5)] 


G 6 [TV(6)J 


2 


(2) 


0.750000 


0.5000 




— 


— 


3 


(3) 


1.310371 


0.6609 








4 


(4) 


1.785826 


0.7715 




— 


— 


5 


(5) 


2.245187 


0.8651 


11 [3] 


— 


— 


6 


(6) 


2.654039 


0.9406 


1.000 [6] 


- 


— 


7 


(7) 


3.064186 


1.0106 


1.000 [5] 


1.000 [2] 


— 


8 


(8) 


3.443409 


1.0714 


0.641 [4] 


0.821 [4] 


- 


9 


(9) 


3.809782 


1.1269 


0.965 [3] 


0.957 [6] 


- 


10 


(10) 


4.164990 


1.1783 


1.000 [2] 


0.861 [8] 


- 




(9,1) 


0.021989 


1.2453 


0.965 [3] 


0.957 [6] 


- 


11 


(11) 


4.513275 


1.2265 


0.940 [2] 


0.894 [8] 


1.000 [1] 




(10,1) 


0.009876 


1.2878 


1.000 [2] 


0.861 [8] 


- 


12 


(12) 


4.838966 


1.2700 


- 


1.000 [12] 


- 




(11,1) 


0.015345 


1.3286 


0.938 [2] 


0.895 [8] 


1.000 [1] 


13 


(12,1) 


5.166798 


1.3659 


- 


1.000 [12] 


- 




(13) 


0.005061 


1.3130 


1.000 [1] 


0.894 [10] 


0.932 [2] 


14 


(13,1) 


5.485915 


1.4033 


1.000 [1] 


0.893 [10] 


0.933 [2] 




(14) 


0.003501 


1.3527 


- 


0.938 [12] 


1.000 [2] 


15 


(HI) 


5.792094 


1.4383 


- 


0.938 [12] 


1.000 [2] 




(15) 


0.009031 


1.3906 


- 


0.885 [12] 


0.963 [3] 


16 


(15,1) 


6.093421 


1.4719 


- 


0.882 [12] 


0.962 [3] 




(16) 


0.012200 


1.4266 


- 


0.897 [12] 


0.993 [4] 




(16) 


0.012635 


1.4267 


- 


0.747 [12] 


0.884 [4] 


17 


(16,1) 


6.388610 


1.5042 


- 


0.891 [12] 


0.993 [4] 




(16,1) 


0.000365 


1.5042 


- 


0.746 [12] 


0.884 [4] 




(17) 


0.015766 


1.4611 


- 


0.738 [12] 


0.810 [5] 


18 


(17,1) 


6.678830 


1.5353 


- 


0.738 [12] 


0.810 [5] 




(18) 


0.018611 


1.4941 


1.000 [2] 


0.829 [8] 


0.920 [8] 


19 


(18 1) 


6.964146 


1.5654 


1.000 [2] 


0.827 [8] 


0.920 [8] 


20 


(19,1) 


7.247181 


1.5946 




0.838 [12] 


0.918 [7] 




(18,2) 


0.004264 


1.6285 


0.991 [2] 


0.824 [8] 


0.913 [8] 


21 


(20,1) 


7.522378 


1.6226 




0.792 [12] 


0.917 [8] 




(19,2) 


0.004668 


1.6557 




0.847 [12] 


0.927 [7] 


22 


(21,1) 


7.795469 


1.6499 


1.000 [1] 


0.877 [10] 


0.880 [10] 




(21,1) 


2.5 ■ 10~ 7 


1.6499 


1.000 [1] 


0.859 [10] 


0.866 [10] 




(20,2) 


0.000976 


1.6821 




0.801 [12] 


0.935 [8] 




(20,2) 


0.001053 


1.6820 




0.763 [12] 


0.909 [8] 



both cases the outermost shell is significantly narrower than the second one and this trend becomes more pronounced 
as TV increases. This is easy to understand because the effect of the confinement is strongest at the cluster boundary, 
i.e. in the outer shell, whereas the inner shells are partially shielded from the trap potential by the surrounding particle 
shells. In contrast, the behavior of the inner shells is not that systematic: in one case (TV = 149) the third shell is of 
similar (relative) width as the second, in the other case (TV = 160) the inner shell is much narrower. The reason are 
symmetry effects which particularly strongly influence the width of the innermost shell (the cluster TV = 160 has a 
closed inner shell with 12 particles which is very narrow). 

In Tab. H]we also provide the first excited states which correspond to metastable shell configurations which are 
different from the ground state. While the overall trend is a rapid decrease of the excitation energy (energy gap to the 
ground state) with increasing TV, some additional systematics is observed. Clusters which open a new shell typically 
possess a close metastable state. For example, for TV = 13 the relative stability of the configurations {TV, 0} and 
{TV — 1, 1} changes, the latter becomes the ground state and the former the first excited state, cf. TabQ] A similar 
trend is observed not only when a new shell is opened but also when an additional particle moves onto the inner shell 
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FIG. 2: Voronoi construction for the cluster N — 17 - the two energetically lowest states with shell configuration N = {1, 16}. 
White (grey) areas are hexagons (pentagons) - indicating the number of nearest neighbors of the corresponding particle (black 
dot), a: ground state, b: first excited ("fine structure") state, c: arrangement of the four particles surrounded by hexagons - 
the two states differ by rotation of the edge AB, black [white] circles correspond to case a) [b)]. 

between the states {Ni — 1, N2} and {iVj., N2 — 1}. Away from these transition points the energy differenc increases 
and eventually another configuration becomes the first excited state. 

An interesting observation is that frequently simulations yielding the same shell configuration resulted in different 
total energies, see e.g. N = 16, 17, 22 in Tab. |T| The differences are much larger than the simulation error, moreover, 
the energies are reproducible. The obvious explanation is that the state of a cluster is not completely determined by 
its shell configuration (as it is the case in 2D). In addition, there exist further excited states, i.e. a "fine structure", 
which are due to a different particle arrangement and symmetry within one shell. To understand the differences in the 
structure of these states with same shell configuration we analyzed the intrashell symmetry by performing a Voronoi 
analysis, i.e. by constructing polygons around a given particle formed by the lines equally bisecting nearest-neighbor 
pairs on the shell, cf. the example of N = 17 shown in Fig. [3 Interestingly, both states do not differ with respect to the 
number of polygons of each kind on the outer shell: there are N(5) — 12 pentagons and N(6) = 4 hexagons. However, 
what is different is the arrangement of the polygons: in one case, the four hexagons form a perfect tetrahedron ABCD 
and are separated from each other by pentagons, cf. Fig. [5Ja, in the other two pairs of hexagons touch, see Fig. |2Jb, 
and the tetrahedron is distorted, as shown in Fig.[2Jc. Two edges remain practically constant (AB w CD w 1.63), but 
the edge AB rotates with respect to the first case by an angle of 34 degrees resulting in a reduction of edges BC and 
AD to about 1.24 while AC and BD increase to 1.94. From this we conclude that of two states the one with the more 
symmetric arrangement of the Voronoi polygons, i.e. (Fig. |2]a), has the lower energy. To quantify this topological 
criterion, we introduce the Voronoi symmetry parameter 



where Nm denotes the number of all particles j in the shell, each of which is surrounded by a Voronoi polygon of order 
M , (M nearest neighbors) and Ojk is the angle between the jth particle and its kth nearest neighbor. A value G5 = 1 
(C?6 = 1) means that all pentagons (hexagons) are perfect, the magnitude of the reduction of Cm below 1 measures 
their distortion. Inspection of the values of Cm for the two {16, 1} configurations for N = 17 (Tab. HJ) reveals that 
the state with lower energy has higher values for both G5 and Gq than the second, confirming our observation above. 
This result is verified for all other N (of course it applies only to states with the same shell configuration) . 

Having obtained with Cm a suitable symmetry measure which is sensitive to the relative stability of ground and 
metastable states, we now return to the issue of the overall cluster stability. To this end we compute the mean Voronoi 
symmetry parameter (MVSP) by averaging over all Gm of a given shell weighted with the respective particle numbers 
N(M). The results for the two outer shells for N < 160 are included in Fig. ^ We clearly see that magic clusters 
have not only a high binding energy but also a prominent symmetry [l^, see in particular N = 12, N = 38, N = 103 



In summary, in this paper we have presented extensive simulation results for spherical Coulomb clusters with 
N < 160. The observed lowest energy states for N > 60 are, in most cases, lower than those previously reported 
and should be valuable reference for experiments with classical Coulomb balls in dusty plasmas or ultracold ions. 
Moreover, the ground state results (shell configurations) are expected to be important also for quantum Coulomb 
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and N = 116. 
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TABLE II: Shell configurations, energy per particle for the lowest lying states, and mean shell radii ri,2,3 [l8| 



N 


Config. 


E/N 




T2 


7*3 


28 


(25 3) 


9.348368 


1.8525 


0.6889 




29 


(25 41 


9 595435 


1.8992 


0.7987 




30 


(26 41 


9 838965 


1.9198 


0.7961 




31 


(27 4) 

' 5^7 


10.079511 


1.9399 


0.7926 




44 


(36 8) 


13.020078 


2.2454 


1.0845 




54 


(44 101 


15 085703 


2.4186 


1.1872 




55 


(43,12) 


15.284703 


2.4618 


1.2772 




56 


(44 121 


15.482144 


2.4743 


1.2770 




57 


(45 121 


15 679350 


2.4869 


1.2763 




58 


(45 12 11 


15.875406 


2.5126 


1.3765 




59 


(46 12 11 


16.070103 


2.5247 


1.3764 




60 


(48 121 


16.263707 


2.5236 


1.2754 




64 


(49 14 11 


17 027289 


2.6101 


1.4478 




65 


(50,14,1) 


17.215361 


2.6212 


1.4477 




80 


(60 19 11 


19 936690 


2 8369 


1.6002 




84 


(61 21 2) 


20 632759 


2.9064 


1.7140 


0.5426 


94 


(67,24,3) 


22.325841 


3.0347 


1.8356 


0.7001 


95 


(67 24 41 


22.491878 


3.0522 


1.8848 


0.8089 


96 


(68 24 41 


22.657271 


3 0606 


1.8846 


8083 


97 


(69,24,4) 


22.822032 


3.0687 


1.8849 


0.8095 


98 


(69,25,4) 


22.986199 


3.0864 


1.9055 


0.8081 


99 


(70,25,4) 


23.149758 


3.0945 


1.9056 


0.8071 


100 


(70,26,4) 


23.312759 


3.1117 


1.9259 


0.8055 


101 


(70 27 41 


23.475164 


3.1291 


1.9450 


0.8028 


103 


(72 27 41 


23 798274 


3.1451 


1.9443 


0.8017 


105 


(73 28 41 


24 120223 


3.1696 


1.9641 


8020 


107 


(75 28 41 


24 439666 


3.1850 


1.9640 


0.8011 


109 


(77,28,4) 


24.757151 


3.2005 


1.9638 


0.8006 


111 


(77,29,5) 


25.072584 


3.2322 


2.0249 


0.8968 


113 


(77,30,6) 


25.385842 


3.2637 


2.0831 


0.9640 


115 


(77 32 61 


25 697308 


3.2949 


2.1162 


9630 


117 


(79 32 61 


26 007089 


3.3094 


2.1158 


0.9622 


119 


(81 32 61 


26.315442 


3.3237 


2.1156 


0.9624 


121 


(83,32,6) 


26.622118 


3.3379 


2.1154 


0.9614 


123 


(83,34,6) 


26.927195 


3.3672 


2.1493 


0.9625 


125 


(84 34 71 


27 230458 


3.3884 


2.1850 


1.0340 


128 


(85 35 81 


27 682123 


3.4235 


2 2358 


1.0922 


130 


(86,36,8) 


27.981234 


3.4445 


2.2501 


1.0917 


133 


(88 37 81 


28.427062 


3.4718 


2.2642 


1.0912 


135 


(88,38,9) 


28.722421 


3.4992 


2.3110 


1.1436 


137 


(90 38 91 


29.016328 


3.5119 


2.3110 


1.1440 


139 


(91 39 91 


29 308774 

*— • / ■ ' >\ ' y..> t t l 


3.5316 


2.3251 


1.1430 


141 


(92,40,9) 


29.599900 


3.5514 


2.3387 


1.1417 


143 


(93,40,10) 


29.889733 


3.5707 


2.3689 


1.1932 


145 


(94 41 101 


30 178106 

' / v7 ■ 1 1 I ' 1 V7 X / 


3 5898 


2.3825 


1.1920 


147 


(95 42 101 


30 465219 


3.6087 


2.3957 


1.1923 


149 


(96,43,10) 


30.750998 


3.6273 


2.4090 


1.1926 


151 


(96,43,12) 


31.035390 


3.6524 


2.4659 


1.2814 


153 


(97,44,12) 


31.318528 


3.6708 


2.4781 


1.2811 


154 


(98,44,12) 


31.459632 


3.6768 


2.4777 


1.2810 


155 


(98,44,12,1) 


31.600488 


3.6887 


2.5042 


1.3846 


156 


(98,45,12,1) 


31.741100 


3.7006 


2.5169 


1.3839 


158 


(100,45,12,1) 


32.021294 


3.7122 


2.5166 


1.3834 


160 


(102,45,12,1) 


32.300405 


3.7238 


2.5161 


1.3833 
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clusters (e.g. in quantum dots) in the strong coupling limit, as for 2D systems it was found that, in most cases, they 
have the same shell configuration as their classical counterpart Further, we have presented an analysis of the 

lowest excited states of small clusters. Besides metastable states with a shell structure different from the ground state 
we identified "fine structure" states which are characterized by different particle arrangement within the shells. These 
states have a lower symmetry which is linked to higher values of the total energy. Despite the decreasing values of the 
excitation energy with increasing N, knowledge of the lowest metastable states is important for understanding the 
dynamic properties of mesoscopic clusters. We expect that the collective excitations of the clusters, i.e. the normal 
modes which are excited in the system if kinetic energy is supplied will be strongly influenced by the metastable 
states. Further, these states sthould be of importance for the melting behavior of mesoscopic Coulomb balls. 
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